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ABSTRACT 


A new inversion method for remote sounding of planetary 
atmospheres is presented which appears to have severed significant 
advantages over the conventional methods. This method is applicable 
to high resolution observations where the spectral lines are f ully 
resolved, and is based on matching the calculated slopes of the 
spectral line profiles with slopes of the observed lineshapes. The method 
is applied to inversion of ozone absorption lines in the earth's 
atmosphere and the results are compared with those obtained by a 
conventional method. The proposed method is seen to provide a 
significant improvement in the overall accuracy of the retrieved 
profiles, with higher vertical resolution and higher levels which 
may be probed. 


1. INTRODUCTION 


Following Kaplan’s suggestion in I959 that atmospheric temperature 
structure and composition distribution may be obtained from an analysis 
of spectral variation in the upwelling thermal radiation, a great deal 
of work has been done in developing this powerful technique. Of the 
available methods today, it is one of the most useful in studying the 
atmospheres of the earth and other planets from ground based observations, 
satellites and space probes. Several review articles and proceedings 
of special conferences covering this subject are now available (e.g. 
Rogers, 1971; Colin, 1972; Deepak, 1977? Twomey, 1977). 

The atmospheric remote sensing technique is based on measuring 
the outcoming radiation at a chosen set of frequencies, 
and finding an inverse solution to the radiative transfer equation. 

This solution may be obtained for vertical temperature profiles as 
well as for the distribution of absorbing gases in the atmosphere. 

In principle, the inverse solution is unique and represents the true 
state of the atmosphere at the time of observations. The accuracy 
of this solution, the vertical resolution and the atmospheric range 
which may be probed, is determined by the observational technique 
(sensitivity, spectral range and resolution) and the inversion method 
employed. Several methods have been proposed and successfully used for 
retrieval of atmospheric profiles from planetary observations (e.g. 
Chahine, 1968, 1970, 1972; Smith, 1970; Conrath, et al. , 1970; 

Conrath, 1972; Prabhakara et al., 197<^ 
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The choice of an inversion method for a given set of observations 
is generally based on the following considerations:^) the sensitivity 
of the solution to a-priori information about the state of the 
atmosphere (ii) the stability and the rate of convergence to the desired 
solution (iii) the effect of system noise on the stability of the 
solution (iv) the number of sampling frequencies used and the spectral 
resolution employed in the observations. 

With high resolution measurements over an appropriately chosen 
spectral range the vertical resolution of the retrieved profiles is 
limited by the width of the weighting functions. When the spectral 
resolution is extremely high (Av « a the half-width) individual 
spectral lines can be completely resolved and the lineshapes determined. 
Atmospheric temperature and composition profiles can then be obtained 
from an analysis of these individual lines. The capability of ultra-high 
resolution sub-doppler spectroscopy is provided by heterodyne 
detection. Well developed in radio astronomy, only recently 
it is finding applications in the infrared [e.g. Teich, 1973; Mumma 
et al., 197^5 1978; Abbas et al. , 1976, 1978 a,b; Betz et al., 1977; 
Menzies and Seals, 1977]. Measurement of stratospheric ozone absorption 
lines with an infrared heterodyne spectrometer and the retrieval of 
vertical concentration profiles have been reported [Abbas et al., 1978a ,b; 
Menzies and Seals, 1978]. The inversion methods used in the analyses 
of the above ozone measurements are based on an inverse solution of 
the radiative transfer equation. 
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The purpose of this paper is to present an inversion method 
for fully resolved spectral lines, which appears to have several 
advantages over the presently used methods . This method assumes 
that the spectral resolution is sufficiently high so that the observed 
line is fully resolved and the lineshape is determined. The 
inversion is then based on finding an inverse solution to the derivative 
of the radiative transfer equation with respect to frequency* This 
implies that the slope of the observed line is matched with the slope 
of the synthetic line computed for the retrieved atmospheric parameters. 
As will be seen, this differential technique leads to (i) narrower 
weighting functions providing a higher vertical resolution (ii) higher 
atmopsheric level which may be probed, (iii) faster and more stable 
convergence and (iv) more accurate retrieved profiles. Thus by taking 
full advantage of the powerful techniques of heterodyne and tunable- 
diode- laser spectroscopy the overall quality of the retrieved atmospheric 
profiles is significantly improved. 

The formulation of the proposed differential method is given 
in section 2. An application of the method to an inversion of synthetic 
absorption line profiles of atmospheric ozone is given in section 3* 

A comparison of some results of this method is made with those of the 
presently used method. 
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2 . Basic Equations 


The emergent spectral intensity 1^ from a non- scattering 
atmosphere is given by the radiative transfer equation which may be 
written in the form 

y t 

X sj (P,T) * B v (T.) t* + J B v (T) K(P,T) dy (1) 

y s 

where B v (T g ) is the Planck function at the source temperature T g , 
s 

t is the transmittance from the surface to the top of the atmosphere, 

St 

the weighting function K(P,T) ~ » and y = --tnP. The first tern 

in Eq, (1) is the surface contribution, whereas the second term represents 
the atmospheric contribution to the spectral intensity. The atmospheric 
transmittance t is 


t v = exp [ - J Et Vi d^], 


( 2 ) 


with k^ as the specific absorption coefficient, and du^ is the element 
of column density for the ith species of the absorbing gas given by 


p T 

du i(P) = q^ (p-) (^r) sec e 2 dz (3) 

o 
v 

where q. is the volume mixing ratio of the gas, 0 is the solar zenith 
X z 

angle, and the subscript o refers to the reference quantities. Equations 
(l)-(3) together vich the equation of state, the hydrostatic equation, 
and a knowledge of the absorption coefficient of the gas (the lineshape 
function and the line strength) provide an analytical relationship 
between the observable quantity 1^ and the atmospheric parameters. 

The spectral intensity I v from a planetary atmosphere is 
measured for an appropriately chosen set of frequencies for which the 
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weighting functions K(F,T) are well distributed over the atmosphere. 

The atmospheric profile may be obtained through an inverse solution 
of equation (1). 

If the temperature profile and the surface pressure is known, 
the concentration profile of the absorbing gas may be inferred. 
Alternatively, if the concentration profile is known, the temperature 
profile and the surface pressure may be obtained. With low resolu- 
tion measurement (Av> a) individual spectral 

lines cannot be resolved and average spectral intensities 
are measured. An inversion of such data for evaluation 
of atmospheric profiles is made over an absorption band, with the 
absorption coefficient, the transmittances and the weighting functions 
gener ally representing an average over several or many lines. 

With ultra-high resolution measurements , such as in heterodyne and 

6 7 

tunable-dlode-laser spectroscopy where resolving powers as high as 10-10* 
may be achieved, the lineshapes of individual spectral lines are 
fully resolved. Assuming that the true lineshape function is known, 
the atmospheric profiles may be obtained through an inversion of one 
or two lines. Inversion of individual spectral lines offers several 
advantages over inversion of low resolution measurements . Since no 
averaging over a number of lines has to be made, the weighting functions 
which may be generated from fully resolved lines are narrower and 
reach higher into the atmosphere. In addition the accuracy of the 
retrieved atmospheric profiles is higher, because the errors due to 
uncertainties in the molecular line parameters have been reduced to 
those for only one or two lines. 
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The theory of lineshapes corresponding to various broadening 
mechanisms is fairly well understood (e.g,, Breen, 1961, Benedict, 

1962; Aller, 1963; Goody, 1964 ). For conditions of planetary atmospheres, 
when no non- thermal processes are assigned (such as maser- laser 
fluorescence emissions and gain narrowing) the two important broadening 
mechanisms in the infrared are pressure or collisional broadening and 
Doppler broadening. The mixed lineshape (Voight lineshape) 
which includes both pressure and 

Doppler broadening effects, is generally representative of lineshapes 
formed in planetary atmospheres. Various modifications to this 
lineshape become necessary only for special considerations, such as 
in contributions from the wings (Benedict’s modification) for frequencies 
(v-v ) >> a, and in high pressure regimes (when the mean free path 
becomes less than the wavelength) where collisional narrowing sets in. 


If the lineshape is completely resolved with ultra-high resolution 
measurements, full advantage of the knowledge of the lineshape function 
may be taken the power of heterodyne method realized f ully by 
employing an alternative method of inversion. This method is based on 
solving for the slope of the spectral intensity. 

The slope is given by the derivative of Eq. (l) with respect to v 


ar * s 

v = W “K + J 

y s 


9v 


b v (t) fjj (p.t) ay 


( 4 ) 


or 


( 5 ) 


K *&> K * I v t) v p * T > * 

jf; : y s 

Where the dot refer*; derivative with respect to v, and it is 

1 v / i!k *Vv ?*'•• '' r ' ' ? 

assumed that the fr equity interval is sufficiently small so that 

dB^/dv » 0. ‘;S* •' 

i’or ground-bastd observations of the earth's atmosphere in the 
solar absorption mode,the atmospheric contribution in (1) and (5) is 
usually negligible, and the observed intensity and its derivative are 


and 


I° b s B (T ) T S 
V V v s' V 


I cb = B (T ) T s 
V v v s' V 


original page is 

OF POOR QUALITY 


( 6 ) 


(7) 


Since Bqs. (5) and (7) are identical in form to the radiative transfer 

a 

equations (1) ar* (6), an inverse solution to 1^ for the temperature or 
gas concentration profiles may be obtained by the usual iterative 
techniques. The differential quantities I v ° b and t® nay be confuted 
from observed spectral line profiles, and are directly measurable 
quantities in systems based on tunable diode lasers. This paper is 
restricted only to inversion of ground-based solar measurements of 
the earths atmosphere for evaluation of gas concentration profiles. 

Temperature sounding of atmospheres based on Eq. (5) will be considered 
in a subsequent paper. 

It should be noted that, in principle, inversion methods based on 
higher order derivations of the radiative transfer equation involving 
d n l/dv n , S n t/bv n etc. could be considered. However, in the absence of any 
new information and the inevitable presence of noise in the observed radiances, 

t 

no significant advantage is likely to be gained. 
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3. Gas Concentration Profile of the Earth's Atmosphere 


The concentration profiles of the earth’s atmosphere may be 
obtained through an inverse solution of Eq* (7)» We have 

P 

T v = expC ‘ i j f k v< ( p ' T) 4 ( p > (”) 

® o * i 


dr , * P dk 

*7 " T v [ - i J o l “ 5 zr < P > T ) 4 i( p > dP ^> 


^ - * * J F K, < P > T > 4 i( p ) «- 


Eq. (8) is similar in form to Eq. (21) of Abbas et al. (1978) 


(8) 


(9) 


T v = * i J o = k Vi dP ( 10 ) 

and may be solved with a similar iterative procedure with k^ as the 
weighting functions. 

The weighting functions k^, in general, are expected to be narrower, 
peaking at higher levels than the corresponding functions k . This may 
be seen in the pressure broadening regime where the absorption coefficient 
is 


* V (P»T) 



S(T) a (P,T) 
( v -v o ) 2 + ^(PjT) 


(H) 


where S(T) is the line strength and a(P,T) is the half-width given by 
a(P,T) = a o (P/P Q ) (■^)' 1 '/ 2 , with the subscript o referring to the reference 
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values* For a constant temperature profile Eq. (H) shows that 

• • • .5 ^ ''ft 

plotted as a function of I oaks at a pressure level where (v-v c ) — 
Oc(P,T) * The weight in^ funct ion k on the other hand is 


ORIGINAL PAGE IS 
Of POOR QUALITY 


L S ( T ) a ( p > T ) (w 0 ) 

vN " C(v-vJ 2 +a a (P,T )f 


( 12 ) 


The form of fiq, (32) indicates that the plots of k^ as a function of 

P are narrower than those for k^. In addition, setting the derivative 

dk^/dos * 0 indicates that in a pressure broadening atmosphere, k^ 

peaks at levels where (v-v 0 ) - a(P,T)/V3« The functions k^ for 

planetary atmospheres, in general, are thus narrower and peak 

at higher levels than the corresponding functions k^. 

« 

The normalized weighting functions k^ and k^ suitable for inversion 
of the atmospheric ozone line (for solar absorption mode) with (J* 

J",K",K. ") identified as (29,3,26;30,3,27) with the line center v = 1011.66 
cm 1 (lower state energy 422.96 cm \ S = 0.4l4 cm’ ^ cm- atm)” 1 , 

Cl = 0.07 cm are shown in Figs. 1-2. A comparison of the two sets 
shows that the plots of k^ are significantly narrower and peak at 
higher levels in the atmosphere. The weighting function k^ for 
v * 1011.670 cm* 1 (Av = O.OOe cm' 1 ), for example, which peaks at 
~ 23 km has a width (FWHM) ~ 18 km. The function on the other 
tend, peaks at -w 28 km, has a width ~ 13 km and is thus expected to 
provide higher vertical resolution. A comparison of the vertical 
resolutions of atmospheric ozone concentration profiles retrieved 
by the two method given in the following section. 
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The iterative solution based on Bq. ( 9 ) is similar 
to that employed in tne literature, (Chahir.e 1972 ; Abbas at al., 1978) 
and is briefly outlined below. The quantity t^/t^ is calculated 
from the observed spectrum for « chosen set of frequencies for which 
the peaks of the weighting functions (J^) are uniformly distributed over 
the atmosphere. For evaluation of concentration profiles it is 
assumed that the temperature profile is known. An initial guess of the 
mixing ratio profile q? ,rn is made (where j denotes the atmospheric level 
and m the iteration) and the quantities t and t^/t are calculated. 

An improved estimate of the mixing ratios at levels corresponding to the 
peaks of the weighting functions is obtained from 


/ * / \ - 

q j ?m+l s J,m 
1 1 (t /t ) x 

v k V 


The complete profile is interpolated from the corrected values and the 

iteration is continued until the calculated valuescf VT / t .. f 1 converge 

k k 

to the observed velues(T /t) ' u over the chosen set of frequencies 

v k 

in a least square nense. The iteration process is stopped when the 


RMS difference in the above two quantities is on the order of the noise 


level of the measur ements. 
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4. Evaluation of the Earth’s Atmospheric Ozone 

We shall now apply the method discussed in the previous section 
for inverting a synthetic absorption line of the earth's atmospheric 
ozone, and make comparions with the conventional method. The 
synthetic spectrum for a midlatitude winter model of the earth's 
atmosphere showing a strong ozone line at v Q - 1011.6670 cm" 1 is 
given in Fig. 3* Following the method outlined in section 3, the 
mixing ratio profile obtained from an inversion of this line using 14 frequencies 
is shown in Fig. 4a. The initial guess profile and model profile 
corresponding to the synthetic absorption line is also shown. The 
retrieved profile has been extrapolated below ~ 7 km and above 
—35 kau A comparison of the retrieved and model profiles indicates 
excellent agreement in the range 7-35 km. The retrieved total ozone 
content is 0.4(32 cm- atm which compares with 0.399 °f the model profile, 
a residual error — 1$. 

For comparison, an inversion of the above ozone line was made 
by using the conventional method (see Abbas et al., 1978) with the 
retrieved profile shewn in Fig. *ft>. The thirteen weighting functions used in this 
case extend from 0-26 km, and the profile has been extrapolated 
below —7 km and above — 26 km. Although the total ozone content 
of the retrieved profile in this case also shows a residual error 
of < 1$ (~ 0.397 cm-atm), the overall accuracy of the retrieved profile 
is much better in the differential method. The number of iterations 
required for convergence in the differential method is ~ 6-8, comparing 
with 15-20 iterations in the conventional method. Any reduction 



In computation time, however , is largely offset by the additional 
time required to calculate the derivatives. 

For a direct comparison of the vertical resolution provided 
by the two methods, a synthetic absorption line generated from an 
assumed ozone profile (Fig. 5a) with a sharp peak (FWHM ~ 9 km) was 
inverted. The initial guess, the assumed model and the retrieved 
profile s using the two methods are shown in Figs. 5. The inverted 
profile in Fig. 5a following the conventional method, approaches the model 
profile, with a somewhat broadened peak and without reproducing the 
dip at ~ 10-15 ton* The calculated total vertical ozone content of 
0.307 cm-atm compares with 0.3004 em-atm of the assumed model profile 
(~£.2$ error). The ozone profile obtained by the differential method 
(Fig* 5b), on the other hand, shows a significantly better overall 
retrieval with almost fully recovered peak of the assumed model 
profile and the dip at — 10-15 ton. lie calculated total 
vertical ozone content of 0.305 cm-atm shows a residual error of ~1.5 $>. 

The proposed method thus clearly provides substantial improvement in the 
accuracy of the retrieved atmospheric profiles and a higher vertical 
resolution. Atmospheric features on the order of a scale height 
(~7 km for the earth) can be fully recovered. The above discussion of the 
vertical resolutions provided by the two methods is of an "empirical" nature. 
More rigorous considerations of the achievable resolution, based on the 
tradeoff analysis of Backus-Gilbert and Conrath (1972) will be given in a 
subsequent paper. 

It should be remarked here that the ozone absorption line of Fig. 3 
analyzed in this paper is by no means the most appropriate, and the retrieved 
ozone profiles of Figs. 4-5 are not the best which may be obtained by 
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either of the two methods* The results presented, here ere for comparison 
purposes only; further improvements a ay be made in the quality of the 
profiles obtained through the two methods. 

An important consideration in the inversion of spectral lines 
Is the effect of noise on the accuracy of retrieved atmospheric 
profiles. In the earth's atmosphere , for species with mixing ratios 
£l0" -lcf 9 (such as CO^O^CH^NO^etc.), high SNR’s (- 100 ) may 
be achieved in the solar absorption mode so that the effect of 
noise is negligible* However, for less abundant species, such as 
ClOjHgOgjCiM’s etc., or for weak lines, the achievable SNR's may be 
small er and the presence of noise may affect the accuracy of the 
calculated atmospheric profiles. 

For inversion of ozone absorption lines considered in this 
paper, the SNR’s are generally quite high and the effect of noise on 
the retrieved profiles is insignificant. In order to evaluate the 
noise effect, however, we may assume a degraded SNR and use a random 
noise generator with Gaussian distribution to add noise to the synthetic 
line profile (Fig. 3). The noisy data may be smoothed by least square fitting 

methods to yield a uniform line profile. Figs. 6a-b show the envelopes of 
the retrieved profiles obtained from inversion of the synthetic line (Fig. 3) 
with superimposed random noise corresponding to an SNR of ~ 10. Whereas the 
effect of noise appears to be comparable in the two retrieval methods for the 
case considered here, the differential method becomes more sensitive 
to noisy data and poor SNR's. The usefulness of the differential 
method then becomes critically dependent on the sophistication of the 
smoothing technique employed in recovering a well structured line profile. 
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The smoothing process, however, results in a degradation of the 
achievable vertical resolution in both cases. 

Conclusions 

The inversion method presented here is applicable to high resolution 
measurements where the spectral lines are fully resolved, and appears to 
have several advantages over the presently used methods . The principal 
feature of the suggested method are (i) a substantially improved vertical 
resolution and an overall higher accuracy with which the absorbing gas and 
temperature profiles may be retrieved (ii) higher vertical levels. which may 
be probed (iii) faster and more stable convergence to the desired solution. 

An application of the method for evaluation of concentration profiles 
of the earth's atmosphere indicates that ozone profiles to heights ~35 km 
may be retrieved with much better accuracy than the conventional methods , 
and vertical features on the order of a scale height may be f ully re- 
produced. The method works well for data with high SNR's (> 10) but 
is more sensitive to noise for lower SNR's. Smoothing techniques which 
yield well structured spectral line profiles may be employed, but the 
advantage of higher vertical resolution tends to be degraded. 

Although the analysis and the examples considered in this paper are 
restricted to evaluation of absorbing gas profiles in the earth’s atmosphere, 
the method is generally applicable to planetary atmospheres for remote 
sounding of both temperature and concentration profiles. Remote sounding 
of atmospheric temperature profiles will be considered in a subsequent 
paper. 
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LIST OF FIGURES 


Fig. 1. Normalized weighting functions K for the ozone line with 
line center v Q = 1011.6670 cm \ for frequencies with 
Av(cm "*■) = (v-v^) given by 

(1) 0.0010 (2) 0.002 (3) 0.003 (4) 0.004 (5) .005 
(6) 0.008 (7) 0.016 (8) 0.018 (9) 0.020 (10) 0.207 
Fig. 2. Normalized weighting functions for the ozone line with 
line center v Q = 1011.6670 cm \ for frequencies with 
AvCcnf 1 ) = (v-v Q ) given by: 

(1) 0.001 (2) 0.002 (3) 0.003 (4) 0.004 (5) 0.005 

(6) 0.008 (7) 0.016 (8) 0.018 (9) 0.032 (10) 0.042 

Fig. 3* Synthetic atmospheric spectrum showing a moderately strong 

absorption line of ozone (29,3,26; 30,3»27), with line 

center v = 1011. 6670 cm \ a = 0.07 cm 1 lower state energy 
00 

E" = 422.96 cm’ 1 . The midlatitude summer model atmosphere 
with one air mass has been assumed. The analyzed spectral 
region is indicated with arrows. 

Fig. 4a. The retrieved volume mixing ratio profile of ozone obtained 
from an inversion of the synthetic line of Fig. 3 with the 
differential method. The initial guess and the nominal model 
profiles are also shown. 

Fig. 4b. The retrieved volume mixing ratio profile of ozone obtained 
from an inversion of the synthetic line of Fig. 3 with the 
conventional method. The initial guess and the nominal model 
profiles are also shown 



Fig. 5a. The retrieved volume mixing ratio profile of ozone obtained 

with the differential method, from an inversion of the synthetic 
line generated from the assumed sharply peaked ozone profile 
shown in this figure. 

Fig 5b. The retrieved volume mixing ratio profile of ozone, 

obtained with the conventional method, from an inversion 
of the synthetic line generated from the assumed 
sharply peaked ozone profile shown in this figure. 

Fig. 6a. The envelope of the retrieved volume mixing ratio profile of 

ozone, obtained with the differential method, from an inversion 
of the synthetic line of Fig. 3 with superimposed random 
noise of Gaussian distribution corresponding to a SNR ~ 10. 

Fig. 6b. The envelope of the retrieved volume mixing ratio profile of 

ozone, obtained with the conventional method, from an inversion 
of the synthetic line of Fig. 3 with superimposed random 
noise of Gaussian distribution corresponding to a SNR ~ 10. 
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